A Factorization Law for Entanglement Decay 
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We present a simple and general factorization law for quantum systems shared by two parties, 
which describes the time evolution of entanglement upon passage of either component through an 
arbitrary noisy channel. The robustness of entanglement-based quantum information processing 
protocols is thus easily and fully characterized by a single quantity. 



Whenever we contemplate the potential technological 
applications of quantum information theory [l[ , from se- 
cure quantum communication over quantum teleporta- 
tion jij to quantum computation we need to worry 
about the unavoidable and detrimental coupling of any 
such quantum device to uncontrolled degrees of free- 
dom - typically lumped together under the label "en- 
vironment". Environment coupling induces decohcr- 
ence [H, 0, 0], i.e., it gradually destroys the phase re- 
lationship between quantum states, and thus their abil- 
ity to interfere. In composite quantum systems, these 
phase relationships (or "coherences" ) are at the origin of 
strong quantum correlations between measurements on 
distinct system constituents - which then are entangled. 
The promises of quantum information technology rely on 
exploring precisely these non-classical correlations. 

Yet, entanglement is not equivalent to many-particle 
coherences: it is an even stronger property, and hard 
to quantify - all commonly accepted entanglement mea- 
sures @, Q are nonlinear functions of the density matrix 
which describes the state of the composite quantum sys- 
tem, and in particular the coherences. While an elab- 
orate theory on the time evolution of quantum states 
under environment coupling is at our hands, virtually 
no general results on entanglement dynamics have been 
stated. Hitherto, the time evolution of entanglement al- 
ways needed to be deduced from the time evolution of the 
state I, [HI El, E2, E, Q El, G3, E3, EH ■ In the present 
Letter, a direct relationship between the initial and final 
entanglement of an arbitrary bipartite state of two qubits 
(basic units of quantum information) subject to incoher- 
ent dynamics in one system component is derived, which, 
as illustrated in Fig. renders the solution of the cor- 
responding state evolution equation obsolete. Our result 
can be directly applied to input/output processes, such as 
gates used in sequential quantum computing. Moreover, 
it allows to infer the evolution of entanglement under cer- 
tain time-continuous influences of the environment, e.g. 
phase- and amplitude damping. 

Let us consider entangled states of qubit pairs, with 
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FIG. 1: From state evolution to entanglement dynamics. 
Hitherto, the time evolution of entanglement C under open 
system dynamics had to be deduced from the time evolution 
of the state \\). However, it turns out that knowlegde of the 
entanglement evolution of the maximally entangled state \4> + ) 
under the channel 1 & $ suffices to establish a direct mapping 
of C(\x)) onto C(p), without detour over the state. 



one qubit being subject to an arbitrary channel $ - which 
may represent the influence of an environment, of a mea- 
surement, or of both. In order to illustrate the situation, 
we consider a source which emits a particle to the left and 
another one to the right. Each particle on its own carries 
one qubit of quantum information (in general a super- 
position or mixture of two basis states |0) and |1)). We 
therefore also refer to the particles as "left" and "right" 
qubit. Let the particles leaving the source be in a pure 
state |x): 



|x) = V^|oo) + VT^J|ii) 



(i) 



with < us < 1, i.e., for values of us between zero and 
one the particles are in a coherent superposition of both 
being in state |0) and both being in state Any pure 
state can be written in this form, modulo local unitary 
operations. Since our results are not affected by these lo- 
cal unitaries, as we will see below, this choice of the pure 
initial state does not impose a restriction of generality. 
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b) 

FIG. 2: Entanglement decay scenarios, a) Laboratory sce- 
nario: one of the qubits of the initial state \x), the "right" 
one, undergoes the action of a general quantum channel $. 
b) Dual scenario: the same end result is obtained by inter- 
changing the role of states and channels: the "left" qubit of 
a mixed state p$ undergoes the action of the quantum chan- 
nel $ x . Circles represent sources of entangled states, squares 
symbolize channels. 




FIG. 3: Quantum Teleportation identity. The teleportation 
protocol transfers the incoming state from the left to the out- 
going state on the right [23| . The procedure consists of a 
composite Bell measurement on the incoming qubit and on 
one qubit of an auxiliary, maximally entangled state \(j> + ). 
We restrict to the case where the measurement results in a 
projection M.^ on the state \4> + }- 



We quantify entanglement by concurrence C [19j |. 
which implies for the state C(|x)) = 2^lu(1 — uj). 
For uj in ([1]) equal to zero or one. the state's entangle- 
ment and hence its concurrence vanishes; uj = 1/2 implies 
|x) = one of the maximally entangled Bell states, 

with maximal concurrence one. 

Now the right qubit traverses an arbitrary quantum 
channel $, as illustrated in Fig. [5^), and we want to de- 
rive the qubits' entanglement hereafter. To do so, note 
that the qubits' final state must be the same as in a dual 
picture |20j, where the roles of initial state and channel 
are interchanged, as depicted in Fig. [2}d). Thus, the two- 
qubit state |x) is identified with a qubit channel $ x , and 
the qubit channel $ with a two-qubit state p$; symboli- 
cally: 

(i®$)ix)(xi = (h®*)p$ _ (2) 

p' p 

Here, p' = Tr[(l ® $) |x)(xl] and p = Tr[($ x ® 1) p$] are 
the probabilities for channels $ and $ x to act on the states 
|x)(xl and p$, respectively. Thus, we also account for 
non-trace-preserving channels, where the particle number 
is not conserved. 

We now need to determine $ x and p$ explicitcly. For 
this purpose, we first remember that Quantum Telepor- 
tation [2j is a means to transfer the state of one system 
to another one, in principle with perfect fidelity. Conse- 
quently, teleporting the right qubit of the state |x) as- 
sisted by the maximally entangled state |</> + ) leaves the 
state |x) invariant. This invariance is depicted in Fig. [3] 
We therefore obtain the same final state as in the situ- 
ation we considered so far (Fig. [2^)) - if we replace the 
source preparing state |x) by a source which prepares |x) 
followed by a teleportation of the right qubit as shown in 
Fig. r?] Let us now consider the source of the qubit pair 
in state \4> + ), which we inserted with the teleportation. 
The succession of processes influencing the left qubit and 



those acting on the right qubit of the pair can be altered 
without consequences for the final state. For this reason 
we can replace the source producing the state \4> + ) to- 
gether with the channel $ acting on the right qubit by 
yet another source which immediately prepares the state 

P$ :=(l®$)|0+)(0+|/p", (3) 

where p" = Trp ® $)|0+)(0+|], see Fig. H The result- 
ing scheme in Fig. [5] transfers entanglement between the 
qubit pairs prepared in states |x) and p$ to entanglement 
between the left qubit of the first pair and the right qubit 
of the second pair. This scheme is called entanglement 
swapping dEH^. 

Finally, we define $ x to be the channel corresponding 
to the change of the left qubit of />$ in Fig. which 
includes a projection of the left qubit of state p$ 

and the right qubit of |x) on \(j> + ). Channel $ x can be 
interpreted as imperfect teleportation assisted by state 
|x), leaves the resulting state in general non- normalized, 
and can be expressed in the particularly simple form: 

($ x ® 1) p $ = (M ® a) p$ (Aft <x> a) , (4) 

with M = (Vw|0)(0| + x/T - U|l)(l|)/V2. The nor- 
malized final state ($ x ® t)p$/p is the same as (a ® 
§)\x){x\/p', as spelled out by © and in Fig. H but 
the entanglement evolution induced by the particular 
channel $ x can be deduced more easily, as we will now 
demonstrate. The concurrence C of the final state 
p' = (a®$)|x)(xl/V is given by 

C(p') = max {O, yfe - - - y/u} , (5) 

where the £i are the eigenvalues of the matrix p' ■ p' , in 
decreasing order, with p' = [a y ® a v ) ■ p'* ■ (a v g) a y ), and 
p'* the complex conjugate of p' , in the canonical basis. 
In order to evaluate this expression, we use relation © 
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FIG. 4: Building in the teleportation identity. The right 
qubit of |x) undergoes the action of a quantum channel $, 
after intermediate teleportation. The maximally entangled 
state \4> + ) together with the action of the channel 1 ® $ yield 
the source of the mixed state p$. 



FIG. 5: Entanglement swapping. Transfer of entanglement 
between the qubit pairs of |x) and p$, respectively, to en- 
tanglement between the outgoing pair of qubits. The Bell 
measurement together with the source of the entangled 

state \x) constitute the quantum channel $ x for the left qubit 
of state p$. 



together with (01 . We write explicitly: 

p' ■ p' = \ (M ® 1) p $ ■ \Mu y M ® <j y ] ■ p* $ ■ [Ma y ® a y ] , 

(6) 

where we employed that M = M+ = M*. For invertable 
M 2J| , it follows that the eigenvalues of p' ■ p' and p$ ■ p~$ 
are proportional, since 



det [p 1 -p'-^t] 



det 
det 



(M ® t) det [M ® 1] det \p' ■ p> - £1 



(M®1) l p'-p'(M® 1) -£1 



(7) 



4p 2 



- lo) det [p$ ■ p$ - pi] 



where p — £ (w(l — o>)/4p 2 ) , and wc used Ma y M = 
\/ — aj)a y /2 in order to obtain the last equality. 
Eq. (J5]), together with the definitions of p$ = (1 ® 

$)\ f)(<t> + \/p ", P' = (If §)\x)(x\/p', and C(\ X )) = 
2y/uj(l — w), thus lead [25( to our central result: 

C[(l®$)|x)<x|] = C[(l®$)|0+)(0+|] C(\ X )) (8) 

- the entanglement reduction under a one-sided noisy 
channel is independent of the initial state |x) and com- 
pletely determined by the channel's action on the maxi- 
mally entangled state. Thus, if we know the time evolu- 
tion of the Bell state's entanglement, we know it for any 
pure initial state [26| . This result can also be interpreted 
in terms of entanglement swapping between a pure state 
|x) and a mixed state p$, leading to the final state p' , due 



to the equivalence of the processes represented in Figs. [2] 
and [5] 

The factorization law |8]) can be generalized for mixed 
initial states po, by virtue of the convexity of entan- 
glement monotones such as concurrence, and given an 
optimal pure state decomposition po = Ej pj \tpj){4>j |> 
in the sense that the average concurrence over this 
pure state decomposition is minimal [27| . It then im- 
mediately follows, by convexity, that C [(1 <g> $) po] = 

C [Ej Pi (1 ® $) Hi) (1>j l] < Ej PjC p ® $) |^-) |], 
and application of ||SJ) leaves us with 

C[(l®$)po]<C7[(l®$)|0+)(0 + |] C(p ). (9) 

This inequality holds for all one-sided channels $ , and has 
an immediate generalization for local two-sided channels 

$!®$ 2 = ($1® i)(i®$ 2 ): 

C[($i®$ 2 )po] <C[($i®l)|0+>(0 + |] (10) 
x C[(l®$ 2 )|0 + )(0+|]C(p o ) • 

The concurrence after passage through a two-sided chan- 
nel is thus bounded from above, which immediately 
implies a sufficient criterion for finite-time disentangle- 
ment n [is, m of arbitrary initial states, in terms of 

the evolution of the concurrence of the maximally en- 
tangled state under either one of the one-sided channels 
(e.g., choose $i or $2 induced by infinite temperature or 
depolarizing environments). 

Let us finally identify relevant cases when equality in 
(IH]) holds. For that purpose, wc consider mixed states 
that are obtained after the application of a one-sided 
channel to an arbitrary pure state, po = (1 ® $) | V'o) (V'o I - 
This occurs, for instance, if the qubit originally pre- 
pared in a pure state suffers amplitude decay, and the 
resulting mixed state again is subject to decay dynamics. 
This is tantamount to the concatenation of channels on 
one side, (1 ® $2) (1 ® $1) \4'o)(4'o\, what can be lumped 
together as one channel which combines both actions, 
(1 ® $ 2 ) (1 ® $1) = l®$2,i- In a similar vein as for (|10p . 
and using the factorization relation |[8]) for pure states, 
wc deduce 

C[(l®$2 4 )|0 + )(0 + |] <C[(l®$ 2 )|0 + )(0+|] (11) 

X C[(1®$!)|0 + ) (</>+|] . 

The initial state's concurrence rescales both sides of the 
equation by the same amount, and therefore it is omit- 
ted in above equation. It is now sufficient to investigate 
the time dependence of the maximally entangled state's 
concurrence under the concatenated channels (much as 
for the evaluation of ([8])): if all of these are of the form 
C(i) = exp(— Ft) (which is the case, e.g., for $1 an am- 
plitude decay and $2 a dephasing channel), then equality 
holds in ©, with $ = % u p = (1 ® $2)^0) <Vo|, and 
C(pt) = cxp(— Tit)C(po) (and, equivalently, for the roles 
of channels 1 and 2 interchanged) . 
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In conclusion, equations (|8I9I10I11|) provide us with the 
first closed expression for the time evolution of a bipartite 
entangled state under general local, single- and two-sided 
channels, without recourse to the time evolution of the 
underlying quantum state itself. This is a general re- 
sult inherited from the Jamiolkowski isomorphism [20| . 
which is here "lifted" from state to entanglement evolu- 
tion (see Fig. [I]), and eases the experimental characteriza- 
tion of entanglement dynamics under unknown channels 
dramatically: instead of exploring the time-dependent 
action of the channel on all initial states, it suffices to 
probe the entanglement evolution of the maximally en- 
tangled state alone. 
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